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Vertical  Directionality  Measurements  of  Ambient  Noise 

in  the  NE  Pacific 

B.  j.  Sotirin  and  W.  S.  Hodgkiss 

Marine  Physical  Laboratory 
Scripps  Institution  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 


I.  Introduction 

The  spatial  distribution  at  low  frequencies  and  low  local  wind  speeds  may  be 
described  as  a  broad  pedestal  of  energy  arriving  between  ±20*  of  horizontal  in  the 
deep  sound  channel  with  levels  up  to  25  dB  above  the  higher  angle  levels  [Ander¬ 
son,  1979;  Burgess  and  Kewley,  1983].  This  distribution  varies  with  local  wind 
speed,  as  the  higher  angle  levels  increase  with  increasing  wind  speed,  undermining 
the  clarity  of  the  pedestal  shape  [Hodgkiss  and  Fisher,  1987;  Sotirin  and  Hodgkiss, 
1989b]. 

The  origin  of  the  sound  arriving  at  the  near  horizontal  angles  warrants  some 
discussion.  Ray  tracing  diagrams  show  that  surface  reflecting  rays  which  do  not 
1  interact  with  the  bottom  arrive  within  a  narrow  band  of  angles  around  ±12*  to  ±15°. 

The  noise  content  of  these  rays  is  enhanced  by  the  sound  generated  by  surface 
ships.  Theoretically,  this  implies  that  the  low  frequency  spatial  distribution  should 
have  a  double-humped  pattern  with  the  width  of  the  notch  between  the  humps 
decreasing  with  depth  [Cavanagh  and  Renner,  1980].  Measurements,  however,  show 
a  pedestal  with  broad  frequency  and  spatial  characteristics.  Two  mechanisms  were 
suggested  by  Anderson  [1979]  to  account  for  the  additional  noise  at  near  horizontal 
angles,  one  based  on  the  shallower  sound  speed  minimum  at  high  latitudes  and  the 
other  on  sloping  bottom  reflection  referred  to  now  as  high  latitude  ducting  and 
down  slope  conversion  (or  the  slope  enhancement  effect)  respectively.  These 
mechanisms  would  convert  noise  traveling  at  near  vertical  angles  to  horizontal 
angles  such  that  the  noise  would  travel  with  low  loss  in  the  deep  sound  channel,  fil¬ 
ling  in  the  predicted  noise  notch.  Levels  within  the  pedestal  of  low  frequency 
ambient  noise  spatial  distribution  then,  are  primarily  from  distant  sources. 

Storm  induced  noise  has  been  observed  at  frequencies  above  150  Hz  as  the 
spectral  levels  increase  with  wind  speed.  At  lower  frequencies  however,  there  has 
been  no  consensus  regarding  the  extent  of  the  wind  dependencies.  Efforts  to  con¬ 
form  to  the  logarithmic  functions  seen  at  high  frequencies  has  led  to  a  variety  of 
predicted  spectral  levels  [Urick.  1984].  Difficulties  arise  in  estimating  the  "wind 
dependent  portion"  [Crouch  and  Burt.  1972J  of  an  omnidirectional  signal  which  may 
be  affected  by  other  sources.  The  obvious  experimental  answer  is  the  use  of  the 
directional  capabilities  of  an  array.  This  has  been  the  approach  of  several  investiga¬ 
tors  [Burgess  and  Kewley.  1983;  Hodgkiss  and  Fisher,  1987]  however  changes  in 
location  and  time  introduce  spectral  level  variations,  which  complicate  the  wind 
induced  noise  estimation. 


The  advent  of  large  aperture  arrays  with  many  elements  provides  the  capabil¬ 
ity  of  measuring  ambient  noise  characteristics  with  high  resolution.  In  this  report, 
use  of  a  large  aperture  array  enables  us  to  present  the  fine-scale  structure  of  the 
ambient  noise  field.  Analysis  and  interpretation  of  this  data  is  documented  in  iSoti- 
rin  and  Hodgkiss,  1989b].  The  array  was  deployed  at  three  nominal  depths  covering 
the  water  column  from  400  m  to  3100  m  at  7.5  m  sampling  increments  from 
Research  Platform  FLIP  which  was  moored  in  the  NE  Pacific.  A  detailed  estimate  of 
the  variability  in  the  spatial  distribution  of  ambient  noise  as  wind  speed  increases 
from  2  m/s  to  12  m/s  over  a  21  hour  period  is  provided  by  the  900  m  array  aper¬ 
ture  which  produces  a  narrow  1*  beam  at  100  Hz. 


IL  Experiment  and  Data  Processing  Description 

The  low  frequency  (20-200  Hz)  vertical  directionality  data  was  collected  by  a 
900  m  aperture  array,  with  120  individual  elements  equally  spaced  at  7.5  m, 
deployed  from  R/P  FLIP.  The  experiment  site  was  about  400  nm  west  of  Monterey 
CA,  at  35*  N,  126°  W  in  about  4700  m  of  water  off  the  Monterey  Fan.  Area  noise  lev¬ 
els  were  relatively  high  [Urick,  1984],  typically  about  80  dB//^Pa/VHz  at  50  Hz.  with 
northern  shipping  lanes  to  the  north,  coastal  shipping  to  the  east,  and  during  the 
experiment,  storm  tracks  to  the  south.  Data,  which  were  recorded  almost  continu¬ 
ously  over  a  20  day  period  during  September  1987,  included  local  wind  velocity  and 
swell  height,  local  shipping  densities,  water  current  velocities  at  43  m  and  90  m 
depths,  array  navigation  data  and  the  low  frequency  acoustic  data. 

The  data  processing  procedure  converted  individual  element  time  series  to 
spatiai  spectra.  The  integrity  of  the  data  selected  were  verified  by  examining  time 
series  and  spectra  for  representative  elements.  A  2  s  time  series  of  the  120  chan¬ 
nels  were  scanned  at  the  beginning,  middle  and  end  of  each  time  series  to  note  any 
channel  variations.  For  each  2  minute  segment  analyzed,  the  time  series  for  a 
representative  element  was  plotted  over  the  segment  to  identify  any  obvious  con¬ 
taminants.  Calculating  the  spectra  of  one  channel  per  array  section  provided  a 
mid-process  verification  of  the  data.  Further  processing  consisted  of  selecting 
specified  frequency  bins,  compensating  for  an  array  sampling  offset  [Sotirin  and 
Hodgkiss.  1989a]  and  beamforming  across  the  array  aperture.  The  data  segment 
(65536  2-ms  samples  for  each  of  the  120  channels)  was  transformed  to  the  fre¬ 
quency  domain  using  a  1024-point  Fast  Fourier  Transform  (FFT)  and  a  Kaiser  Bessel 
(a-2.5)  window.  Estimates  (0.49  Hz/bin)  were  selected  at  5  Hz  increments  from  15 
to  130  Hz,  avoiding  narrow  band  contaminants  at  60,  115  and  120  Hz  by  shifting 
frequency  slightly  to  61.0,  113.8  and  118.7  Hz.  A  frequency  dependent  phase 
correction  was  implemented  to  account  for  the  sampling  offset  and  the  spatial 
conversion  was  accomplished  using  an  FFF  beamformer.  The  beamfonner  was 
implemented  with  M-512  points,  a  Kaiser-Bessel  (a- 1.5)  shading  function  w(n)  and 
a  spatial  sampling  frequency  of  0.131579  samples/meter.  This  sampling  deviates 
slightly  from  the  array  design  element  spacing  of  7.5  m  due  to  strain  in  the  kevlar 
strength  member  under  tension.  The  spatial  samples  were  incoherently  averaged, 
and  converted  from  electrical  angle  <j  to  physical  angle  6  =  asin (<i>/kd).  The  averaged 
magnitude  square  results  were  calibrated  to  dB  re  1  /iPa/VHz  Deg  by  normalizing 

1  ,  . 

by  the  number  of  points  in  the  FFT,  the  power  in  the  window  function  ( —  Y]  w  («)) 

N  __ 

and  the  spatial  sampling  frequency,  and  by  converting  wavenumber  to  degrees  with 
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the  factor  dk^/ds  — - cos( - 9),  where  is  the  digital  spatial  frequency  —  sin0 

180c  180  X 

analogous  to  the  digital  time  frequency  f.  This  normalization  presents  the  data  in 
terms  of  average  narrowband  Intensity  received  by  the  array  per  vertical  degree. 


in.  Results 

The  ambient  noise  directionality  data  presented  were  selected  from  several 
areas  of  interest.  The  passage  of  a  storm  increased  local  wind  speed  from  2  to  12 
m/s  over  a  21  hour  period  from  Julian  day  (Jday)  257  at  2300  Greenwich  Mean  Time 
(GMT)  until  1940  GMT  the  following  day  (Jday  258).  Two  minute  samples  were 
selected  every  20  minutes  during  this  event  to  provide  a  detailed  investigation  of 
the  directional  variability  with  increasing  wind  speed.  A  second  time  series  is 
presented  to  show  the  variability  in  the  spatial  distribution  on  a  shorter  time  scale 
(about  1  hour).  During  a  second  storm,  the  array  was  deployed  at  three  nominal 
depths  which  provides  an  opportunity  to  examine  the  variation  of  wind  effects  with 
depth  at  three  nominal  wind  speeds. 

The  long  time  series  (21  hours)  is  presented  first  as  omnidirectional  levels  as  a 
function  of  frequency  (1/2  Hz  bands  centered  at  15,  20,  25.  30,  35,  40,  45,  50,  55,  60, 
65,  70,  75.  80,  85.  90.  95.  100,  105.  110.  115,  120,  125  and  130  Hz)  and  time  (sampled 
every  20  minutes).  Selected  array  channel  estimates  were  averaged  and  the  mean 
was  removed  from  each  of  the  narrowband  time  series  which  were  then  passed 
through  a  4-sample  (time  corrected)  running  mean  filter  to  provide  the  estimates  in 
Figure  1.  The  corresponding  filtered  wind  speed  estimates  are  also  shown.  This 
same  data  is  presented  in  Figure  2  as  spatial  spectra  versus  frequency,  at  each  20- 
minute  interval  and  corresponding  local  wind  speed,  and  in  Figure  3  as  spatial  spec¬ 
tra  versus  time,  for  equally  spaced  20-minute  samples,  and  for  the  same  24  frequen¬ 
cies.  Although  the  spectra,  calibrated  to  dB  re  1  /<Pa/  per  VHz  per  vertical  degree 
with  positive  angles  looking  up  and  negative  angles  looking  down,  are  plotted  as  a 
function  of  time  during  the  storm,  the  higher  angle  levels  reflect  the  increasing 
wind  speeds.  To  observe  this  variation  with  wind  speed,  selected  angles  (2,  8,  14, 
20,  ±50.  60,  70  and  ±80")  were  extracted,  normalized  and  plotted  individually  in  Fig¬ 
ure  4  at  15,  35,  55,  75,  and  95  Hz.  The  number  to  the  right  of  each  directional  time 
series  is  the  normalization  factor.  The  lower  time  series  represents  the  correspond¬ 
ing  unfiitered  wind  speed  estimates  in  meters  per  second. 

The  second  set  of  time  series  consists  of  consecutive  2-minute  estimates  of 
the  spatial  distribution  to  show  the  short-time  scale  (about  1  hour)  temporal  varia¬ 
bility  at  specific  frequencies  (15,  35,  55,  75,  95  and  1 15  Hz).  Data  acquired  on  Julian 
day  268  at  1000  GMT  with  the  array  deployed  at  a  mid-depth  of  2650  m  and  a  local 
wind  speed  of  2-3  m/s  is  presented  in  Figure  5.  Data  acquired  on  Julian  day  258  at 
1430  GMT  with  the  array  deployed  at  a  mid-depth  of  850  m  and  a  local  wind  speed 
of  10-1 1  m/s  is  presented  in  Figure  6. 

The  variation  in  spatial  distribution  as  a  function  of  depth  is  presented  in  Fig¬ 
ures  7-9.  Nine  2  minute  data  segments  were  processed  with  characteristics  shown 
in  Table  1.  The  mid-array  depth  for  the  three  deployments  was  850,  1750  and  2650 
m.  The  recorded  wind  speeds  corresponding  to  the  data  segments  were  within  0.5 
m/s  of  3.1,  7.2  and  11.3  m/s,  for  the  low,  medium  and  high  designations  respec¬ 
tively.  The  spatial  distribution  is  plotted  at  each  of  24  frequencies  (15,  20,  25,  30, 
35.  40.  45.  50,  55,  60,  65.  70,  75,  80.  85,  90,  95,  100,  105,  110,  115,  120,  125  and  130 
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Hz)  as  a  function  of  array  depth  and  wind  speed.  Figure  7  is  plotted  with  the  data 
corresponding  to  the  shallow  array  depth  in  the  foreground,  followed  by  that  for 
the  mid  array  depth  and  deep  array  depth  each  offset  by  5  mm.  Figure  8  is  plotted 
with  the  data  corresponding  to  the  low  wind  speed  in  the  foreground,  followed  by 
that  for  the  mid  wind  speed  and  high  wind  speed  each  offset  by  5  mm.  Figure  9  is 
the  spatial  distribution  as  a  function  of  frequency  at  a  specific  array  depth  and 
wind  speed. 
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Figure  1 


Omnidirectional  narrowband  21  hour  time  series. 


Figure  2. 

Spatial  spectra  versus  frequency  at  specified  times. 
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Figure  4. 


Directional  2 1  hour  time  series  at  specified  frequencies. 
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Figure  6. 


Spatial  spectra  1  hour  time  series  at  high  wind  speeds. 
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Figure  7. 


Spatial  spectra  versus  array  depth  at  specified  frequencies  and  wind  speeds. 
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Spatial  spectra  versus  wind  speed  at  specified  frequencies  and  array  depths. 
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Figure  9 


Spatial  spectra  versus  frequency  at  specified  array  depths  and  wind  speeds. 
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Figure  9. 


Spatial  spectra  versus  frequency  at  specified  array  depths  and  wind  speeds. 
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